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Abstract 

We present preliminary results of the analyses of B + — > h + ir° (with h + = ir + ,K + ) and B° — » K°tt° 
decays from a sample of approximately 88 million BB pairs collected by the BABAR detector at the 
PEP-II asymmetric-energy B Factory at SLAC. We measure the 7r + 7r° branching fraction and we 
obtain 

B{B + -» vr+vr ) = (5.5^;° ± 0.6) x 10~ 6 

with a significance of 7.7a including systematic uncertainties. We measure the K + ir° and K°ir° 
branching fractions to be B(B + K + ir°) = (12. 8~t\^ ± 1.0) x 10~ 6 and B(B° -* K°ir°) = 
(10.4 ± 1.5 ± 0.8) x 10 -6 . At the same time, the direct CP-violating asymmetries are investigated 
and we find A n +^o = — 0.03±g;Jf ±0.02, A K+n o = -0.09±0.09±0.01 and Akw = 0.03±0.36±0.09, 
where the errors are statistical and systematic, respectively. 
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1 Introduction 



The study of B meson decays into charmless hadronic final states plays an important role in the 
understanding of CP violation in the B system. Measurements of the CP-violating asymmetry in 
the ir + ir~ decay mode can provide information on the angle a of the Unitarity Triangle. However, 
in contrast to the theoretically clean determination of the angle (3 in B decays to charmonium final 
states Jl],|2| the extraction of a from 7r + 7r~ decays is complicated by the interference of b — > uW~ 
tree and b — ► dg penguin amplitudes. Since these amplitudes have similar magnitude but carry 
different weak phases, additional measurements of the isospin-related decays^, B + — > 7r + 7r° and 
B° — ► 7r°7r°, are required to provide a way of measuring the angle a [||. The measurement of the 
branching fraction of the B + — ► 7r + 7r° decay is, in fact, a crucial ingredient, since it is a pure tree 
decay to a very good approximation. Therefore, in this channel direct CP violation, detected as a 
charge asymmetry ( w 4 7r + 7r o), is expected to be zero. 

On the other hand, measurements of B — > Ktt decays can be related to a model dependent 
extraction of the weak phase 7 with a global fit to the observables. In order to do this, several 
models have been proposed, based on different dynamical assumptions for B decays to two light 
pseudoscalar mesons 0-|7f|. All these approaches are able to reproduce experimental values for 
branching fractions, but they do not show a good sensitivity to the value of the weak phase 7. 
Providing information on CP-violating asymmetries could contribute to increase this sensitivity, 
clarifying the theoretical framework and improving our ability to constrain the Unitarity Triangle 
in the (p, fj) plane ||. 

In the case of K°tt°, the ideal measurement is a time dependent analysis of the K^t: final state, 
which is a CP eigenstate. In this way, one would be sensitive to the direct CP-violating asymmetry 
(related to the coefficient of the cosine term) and to the weak phase entering the decay (related to 
the coefficient of the sine term) Q. In the present analysis, we perform a time integrated study, 
providing the direct CP- violating asymmetry defined as 

\A\ 2 -\A\ 2 

A C p 



\A\ 2 + \A\ 2 



where A (A) is the decay amplitude (its CP conjugated) taken into account. We actually investigate 
A K o T o which, in the SM, is equal to .A^o^o, neglecting contributions with more than one weak 



1 s 
boson 



We present here results on the B + -> 7r + 7r°, B + -> K + ir° and B° -> K°tt° decays. The BABAR 
collaboration has previously published (T^j observations of these channels: we now have reduced 
the errors on branching fractions and investigated direct CP-violating effects. 



2 The BAR4R Detector and Data Set 

The data used in these analyses were collected with the BABAR detector at the PEP-II e + e~ 
storage ring. The sample corresponds to an integrated luminosity of about 81 fb _1 accumulated 
on the T(4S) resonance ( "on-resonance" ) and about 9fb _1 accumulated at a center-of-mass (CM) 
energy about 40MeV below the T(4S) resonance ("off-resonance"), which are used for background 
studies. The on-resonance sample corresponds to (87.9 ± 1.0) x 10 6 BB pairs. The collider is 
operated with asymmetric beam energies, producing a boost (^7 = 0.55) of the T(4S) along the 
collision axis. 

1 Charge conjugate modes are assumed throughout this paper unless explicitly stated. 
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BABAR is a solenoidal detector optimized for the asymmetric beam configuration at PEP-II and 
is described in detail in Ref. [jll). Charged particle (track) momenta are measured in a tracking 
system consisting of a 5-layer, double-sided, silicon vertex tracker and a 40-layer drift chamber 
filled with a gas mixture of helium and isobutane, both operating within a 1.5 T superconducting 
solenoidal magnet. Photon candidates are selected as local maxima of deposited energy in an elec- 
tromagnetic calorimeter (EMC) consisting of 6580 CsI(Tl) crystals arranged in barrel and forward 
endcap subdetectors. In this analysis, tracks are identified as pions or kaons by the Cherenkov 
angle 8 C measured by a detector of internally reflected Cherenkov light (DIRC). The DIRC system 
is a unique type of Cherenkov detector that relies on total internal reflection within the radiating 
volumes (quartz bars) to deliver the Cherenkov light outside the tracking and magnetic volumes, 
where the Cherenkov ring is imaged by an array of ~ 11000 photomultiplier tubes. The Cherenkov 
angle 9 C is measured with a typical resolution of 3 mrad, with a separation between kaons and 
pions of 8a (2.5a) at 2GeV/c (4GeV/c). Good separation at high momenta is essential for two- 
body B decay since the boost increases the momentum range of the decay products from a narrow 
distribution centered near 2.6GeV/c in the CM to a broad distribution extending from 1.7 up to 
4.3GeV/c. 



3 Event Selection, 7r° and K° Reconstruction 

Hadronic events are selected based on track multiplicity and event topology. Backgrounds from non- 



hadronic events are reduced by requiring the ratio of the second to zeroth Fox- Wolfram moment [12| 
to be less than 0.95 and the sphericity [ff3| of the event to be greater than 0.01. 

Charged ir and K candidates (except for K® daughters) are reconstructed within the tracking 
fiducial volume and quality criteria are imposed: they are required to originate within 1.5 cm in 
the xy plane and 10 cm in z from the interaction point, to have at least 12 measured drift chamber 
hits and a minimum transverse momentum of 0.1 GeV/c. 

Candidate ir° mesons are reconstructed as pairs of photons with an invariant mass within 3a 
of the nominal ir° mass [14|, where the resolution a is about 8 MeV/c 2 . Photon candidates are 



selected as showers in the EMC that have the expected lateral shape, are not matched to a charged 
track, and have a minimum energy of 30 MeV. The tt° candidates are then kinematically fitted 
with their mass constrained to the 7r° nominal mass. 

mesons are detected in the mode K° — > — > -k + tt . K® candidates are reconstructed 
from pairs of oppositely charged tracks that form a well-measured vertex and have an invariant 
mass within 11.2 MeV/ c 2 (which corresponds to 3.5<r) of the nominal K® mass [14|. The measured 



proper decay time of the K® candidate is required to exceed 5 times its uncertainty. 



4 B Reconstruction 

Charged (neutral) B meson candidates are reconstructed by combining a 7r° candidate with a track 
h + (a Kg candidate). The kinematic constraints provided by the T(45) initial state and knowledge 
of the beam energies are exploited to efficiently identify B candidates. We define a beam-energy 

substituted mass mES = \j Pb ' wnere E\> = (s/^+Pi'Pb)/^; y/s and E\ are the total energies 
of the e + e~ system in the CM and lab frames, respectively, and pi and pb are the momentum 
vectors in the lab frame of the e + e~ system and the B candidate, respectively. An additional 
kinematic parameter AE is defined as the difference between the energy of the B candidate and 
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half the energy of the e + e~ system, computed in the CM system. The toes resolution is dominated 
by the beam energy spread, while for AE the main contribution comes from the measurement 
of particle energies in the detector. These two variables are therefore substantially uncorrelated 
(< 6%). 

For all the decay modes, the signal energy-substituted mass is parameterized on simulated 
signal events by a Crystal Ballf] function. The resolution is found to be about 2.9 MeV/c 2 and it is 
validated by comparing data and Monte Carlo resolutions for decays into open charm final states 
with large branching fractions, such as B~~ — > D°p~, (with p~ — ► 7r~7r and D° — > K~ir + ). 

The AE distribution for the signal /i + 7r° and K°7r° events is described by another Crystal Ball 
function. The mean value of this distribution is directly obtained from the fit for the h + Tr° sample 
and assumed to be the same for K°ir° since the dominant effect comes from the tt° energy scale. 
For K + ir° candidates, the mean of AE is also shifted because the pion mass is assumed for all 
charged tracks in order to extract the yields of both modes and the AE mean from one fit. This 
shifted mean value can be expressed as 

(AE) = - 7boost x (Jm 2 k +p* - y/M*+p^ , 

where p is the momentum of the assumed kaon track, and M„ and Mk are the pion and kaon mass 
values, respectively. We estimate the resolution on AE to be about 42 MeV, based on simulated 
B + -> h + ir° and B° -» K°tt° events and cross-checked on the B — ► D°p control sample. 

Candidates are selected in the range 5.2 < tties < 5.3GeV/c 2 . Different requirements on AE 
specific to each analysis are described later. 

5 Background Rejection 

The dominant background to these channels is the continuum background, coming from random 
combinations of a true 7r° with a track (a true produced in e + e~ — ► qq continuum events 

(where q = u, d, s or c). Another source of background originates from B decays into other light 
charmless meson final states (charmless background). The main contribution to this background 
comes from B + — ► p + it° and B° — ► p^n^ in the case of B + — > 7r + 7r° and B + — ► K* + ir° in the case of 
B+ K + tt° and B° -> K°tt°. Detailed Monte Carlo simulation, off- resonance, and on-resonance 
data are used to study backgrounds. 

In the CM frame the continuum background typically exhibits a two-jet structure, in contrast 
to the isotropic decay of BB pairs produced in T(45) decays. We exploit the topology difference 
between signal and background by making use of two event-shape quantities. The first variable 
is the angle 9 S between the sphericity axes of the B candidate and of the remaining tracks and 
photons in the event. The distribution of | cos# s | in the CM frame is strongly peaked near 1 for 
continuum events and is approximately uniform for BB events. We require |cos# s | < 0.8. The 
second quantity is a Fisher discriminant constructed from the quantities Lq and L2, where Lj 
is: 

Lj = y^Pi\ cos 9i\ j 

i 

and pi (Oi) is the momentum (the angle with respect to the thrust axis of the B candidate in the CM 
frame) of each charged track and neutral cluster not used to reconstruct the candidate B meson. 

2 The Crystal Ball function pa] is a core Gaussian with a power law to describe the left tail. 
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Monte Carlo samples are used to obtain the values of the Fisher coefficients, which are determined 
by maximizing the statistical separation between signal and background events. No requirement is 
applied on T\ instead the distributions for signal and background events are included in a maximum 
likelihood fit as described in the next section. 

On the other hand, charmless background events tend to peak in tties, as do signal events, 
and have the same T distribution as signal, since they are true B decays. Nevertheless they are 
characterized by lower AE values, since at least a pion is lost in the B reconstruction. A cut on 
AE is the only way to reduce this background to a negligible level. We use on-resonance data in 
the negative AE sideband region (—0.40 < AE < — 0.20GeV) to estimate the magnitude of this 
background. The efficiency determined from Monte Carlo events is used to scale the number of 
events in the sidebands and get the expected number of events in the chosen signal region. We 
finally require -0.11 < AE < 0.15 GeV for /i+vr and -0.15 < AE < 0.15 GeV for K°ir°, reducing 
the charmless background contribution to the level of less than 1%. 

A total of 21752 candidates in the on-resonance data satisfy the h + ir° selection, and a total of 
2668 candidates satisfy the K®7r° selection. These two samples enter into two separate maximum 
likelihood fits. 

The final signal selection efficiency e is (26.1±1.7)% for B + -► vr+vr and (21.5±1.5)% for B + -► 
K+ir° events, while it is (28.0 ± 2.0)% for B° -► K°ir° events. The errors on the efficiencies are 
statistical and systematic, combined in quadrature. The dominant contribution to the systematic 
error is due to the imperfect knowledge of 7r° and Kg reconstruction efficiencies (5% and 3% relative 
errors, respectively). The hierarchy of efficiency values comes from the difference in AE lower cut 
(/i + 7r° vs. KgTT ) and from the use of the tt mass hypothesis in the AE calculation (7r + 7r° vs. 

K+TT°). 



6 Signal Extraction 

For each topology (h + 7r° and i^°7r°), an unbinned maximum likelihood fit determines the signal 
and background yields rii (i = 1 to M, where M is the total number of signal and background 
types) and CP asymmetries. The measured asymmetry is defined as: 



A 1 = 



m + ni 



where fii is the fitted number of i th type h~ir° [K°7r°] events and corresponds to /i + 7r° [K°7r°] 
events. The input variables (xj) to the fit are toeS; AE, T and, in the case of charged modes, the 
Cherenkov angle 6 C of the track from the candidate B decay to distinguish between the final states 
with h = tt and h = K. 

The h + TT° extended likelihood function C is defined as 



M 



N 



cxp 



e* n 



M i 



where qj is the charge of the track h in the j th event and, in this case, M is equal to 4 including 
signal and background 7r + 7r° and K + ir°. The M probabilities Vi(xj;cti) are evaluated as the 
product of probability density functions (PDFs) for each of the independent variables Xj, given 
the set of parameters Si which define the PDF shapes. Monte Carlo simulation is used to validate 
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the assumption that the fit variables are uncor related. The exponential factor in the likelihood 
accounts for Poisson fluctuations in the total number of observed events N. 

For the K°tt° mode we need to measure the flavor of the B candidate in order to extract the CP 
asymmetry. We use B flavor tagging information to distinguish B° from B° decays. B candidates 
are defined B° when the other B is recognized to be a B° and vice versa. Details of BABAR B 
flavor tagging can be found in Ref. ||. We have four different tagging categories (k = 1, .., 4), with 
different tagging efficiencies e k and wrong tag fractions w k . We also use untagged events (k = 0) 
in the fit. The extended likelihood is 



/ M \ N 
C = exp I - ^2 n i ) II 

\ j=l / j=l 



AM 



2 

,fe=0 i=l 



where Cj and fj are the tagging category and the measured flavor of the j th event, n k = e k rii is the 
number of events of the i th type in the k th category and M, in this case, is equal to 2 including 
signal and background. The V k (xj-,ai k ) are in principle category dependent. We found that such 
dependence can be ignored without any bias in the fit, for all the PDFs, except for background 
mES and background T. 



The ?Ti Eg shape in background is parameterized by the ARGUS threshold function [16 



f(x) oc xyl — x 2 exp[— £(1 — x 2 )] 

where x = mEs/ m o an d ttiq is the average CM beam energy, and £ is the parameter determining 
the shape of the distribution and is left free to float in the likelihood fits. The background shape 
in AE is parameterized as a second-order polynomial whose coefficients are taken from a fit to the 
on-resonance ttt-es sideband region. The signal distributions have been already described in Sect. |||. 

Events from simulated signal decays and from on-resonance 7ti.es sideband regions are used 
to parameterize the Fisher discriminant PDFs for signal and background events as a bifurcated 
Gaussian^ and a sum of two Gaussians, respectively. Alternative parameterizations for J-, obtained 
from off-resonance data (for background) and from a large sample of fully reconstructed B° — ► 
nir decays (for signal), are used to estimate systematic uncertainties. The 9 C PDFs are derived 
from kaon and pion tracks in the momentum range of interest from a sample of D* + — > D°7r + 
(D° — ► K~tt + ) decays. This control sample is used to parameterize the 6 C resolution (ag c ) as a 
function of track polar angle; double- Gaussian PDFs are constructed from the difference between 
measured and expected values of 9 C for the pion or kaon hypothesis, normalized by the error ag c . 
Tagging efficiencies and mistag fractions are estimated from a sample of fully reconstructed B 
decays. 

The results of the fit are summarized in the first column of Table [l], where the statistical error 
for each mode corresponds to a 68% confidence level interval and is given by the change in signal 
yield raj that corresponds to a —2\nC increase of one unit. In the case of the 7r + 7r° final state, we 
evaluate how the imperfect knowledge of the PDF shapes can affect the significance of the signal. 
We recalculate the square root of the change in — 21n£ with n^+^o fixed to zero for the worst case 
PDF variations and we find a statistical significance of 7.7a for the signal. 

In order to increase the relative fraction of signal events of a given type for display purpose 
we choose events passing requirements on probability ratios. These probability ratios are defined 

3 The bifurcated Gaussian is an asymmetric Gaussian having two different a, one for x > /i and another for x < /i, 
where \x is the mean value. 
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as 7£ s ig = ^gVs/ Yli'Pi arm T^k = "Phi Yls'Psi where ^ s denotes the sum over the probabilities 
for signal hypotheses only, denotes the sum over all the probabilities (signal and background), 
and Vk denotes the probability for signal hypothesis K + ir° (for /i + 7r° only). These probabilities 
are constructed from all the PDFs except that describing the plotted variable. Figures [l] show 
the distributions in ttt-es and AE for events passing all such selection criteria. The likelihood fit 
projections, scaled by the relative efficiencies for the probability ratio requirements, are overlaid on 
each distribution. Since the sample projections in 7ties an d AE are obtained with requirements on 
different probability ratios, the number of signal events appearing in the two projections are not 
the same. The efficiencies for the t«es plots are: 24%, 53% and 65% for 7r + 7r°, K + ir° and K°tt° 
signal events, respectively. For h + ir° AE projection plots, we show a wider window (—0.200 < 
AE < 0.150 GeV) with respect to the signal region used in the maximum likelihood fit (—0.110 < 
AE < 0.150 GeV) in order to show the BB background present at low AE values. The efficiencies 
for the AE plots are: 35%, 48% and 98% for 7r + 7r°, K + tt° and K°7r° signal events, respectively. 

7 Branching Fraction and Direct CP Asymmetry Results 

The branching fractions are defined as 

1 n h+7T o 



B(h + ir°) 



B(ir° -^77) e h+n o -N BS ' 



£ (vr o ^ 77) . B{K o ^ K s) . B{K o ^ . €k>0 . Nb _ > 

where n h +^o (n^o^o) is the signal yield from the fit and e^+^o (e^-o^o) is the reconstruction efficiency 
for the mode h + n° (K°ir°) in the detected vr° (K°) decay chain. N BB = (87.9 ± 1.0) x 10 6 is the 
total number of BB pairs in our dataset. B(-k° — ► 77), B(K° — ► Kg), and B(K® — > tt + tt~) are 
taken to be equal to 0.9880, 0.5 and 0.686, respectively [Q. Implicit in the above equations is the 
assumption of equal branching fractions for T(4S) — > B°B° and T(45) — > B + B~ . 
The CP asymmetry in the K°7r° channel is defined as: 

<Ai = ■A m • (1 + x d ) 

where x d = Am d /T = 0.755 ±0.015 0. Since the flavor of the signal B cannot be determined, we 
apply the correction factor (1 + x^) to take into account the B° — B° mixing and to translate the 
measured asymmetry into the direct CP asymmetry Acp ignoring CP-violating effects in mixing. 
In the case of the charged B, A h+n o = A^ w ° . Table | summarizes the results on the branching 
fractions and the CP asymmetries: the 90% confidence level (CL) intervals for the asymmetries are 
also given. 

Systematic uncertainties on the branching fractions arise primarily from uncertainty on the 
final selection efficiency and uncertainty on due to imperfect knowledge of the PDF shapes. The 
latter is estimated either by varying the PDF parameters within la of their measured uncertainties 
or by substituting alternative PDFs from independent control samples. For K°ir° analysis, mistag 
fractions and tagging efficiencies are varied by la of their measured uncertainties. 

In the /i + 7r° analysis the most relevant systematic uncertainties on the signal yields are due 
to the Fisher shape for both signal and background events (about 5% each), while for the K°tt° 
analysis the AE background and Fisher signal shapes contribute with the largest errors. We 
estimate the systematic uncertainty on the signal yields due to the residual presence of charmless 
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Figure 1: Distributions of tbes (left) and AE (right) for 7r + 7r° events (top), K + ir events (center) 
and K°ir° events (bottom), after additional requirements on probability ratios, based on all vari- 
ables except the one being plotted. Solid curves represent projections of the complete maximum 
likelihood fit result; dotted curves represent the background contribution. For the ities vr + 7r° plot, 

plots, hatched areas rep- 



K + tt° cross-feed is also shown with a dashed curve. For AE ir' 

o 



-7T° (K+7T° 



resent continuum background only. For the AE tt + 7t° plot, 7r + 7r° signal, K + tt° cross-feed and BB 
background are represented by the dot-dashed, dotted and dashed curves, respectively. The BB 
background is clearly present at low AE values. The vertical line represents the AE requirement 
applied in the analysis. 
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Table 1: Summary of fitted signal yields, measured branching fraction B and CP asymmetries Ai 
The first error is statistical and the second is systematic. 



Mode 


Signal Yield 


6(10- e ) 




Ai (90% CL) 


vr+vr 


1251^ ± 10 


5.5+^ ±0.6 


-0.03^;^ ± 0.02 


[-0.32,0.27] 


K+7T 


2391^ ± 6 


12.8±i;f ± 1.0 


-0.09 ±0.09 ±0.01 


[-0.24,0.06] 


K°TT° 


86 ± 13 ±3 


10.4 ± 1.5 ±0.8 


0.03 ± 0.36 ± 0.09 


[-0.58,0.64] 



backgrounds with toy Monte Carlo techniques and we find that it is negligible compared with the 
other effects. 

Systematic uncertainties on the CP asymmetries are evaluated from PDF variations. For 
charged modes, this contribution is added in quadrature with the limit on intrinsic charge bias 
in the detector (0.01). For the K°7r° mode, an additional contribution (0.02) coming from tagging 
has been evaluated by moving the tagging parameters (efficiencies and mistag fractions) by la from 
their nominal value. 

In conclusion, we have presented preliminary measurements of the branching fractions of B + — ► 
7r + 7r , B + — ► K + ir° and B° — > K°ir°. These measurements supersede our previous results. We do 
not observe any evidence of direct CP asymmetry in these channels. 
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